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^\ ' We propose a lattice model, in both one- and multidimensional versions, which may give rise to 

matching conditions necessary for the generation of solitons through the second-harmonic generation. 
The model describes an array of linearly coupled two-component dipoles in an anisotropic nonlinear 
host medium. Unlike this discrete system, its continuum counterpart gives rise to the matching 

■ conditions only in a trivial degenerate situation. A system of nonlinear evolution equations for 
' slowly varying envelope functions of the resonantly coupled fundamental- and second-harmonic 

" waves is derived. In the one-dimensional case, it coincides with the standard system known in 

nonlinear optics, which gives rise to stable solitons. In the multidimensional case, the system prove 
' to be more general than its counterpart in optics, because of the anisotropy of the underlying lattice 

, model. 
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I. INTRODUCTION AND FORMULATION OF 
THE MODEL 



Solitons supported by quadratic nonlinearities, fre- 
quently called "quadratic solitons" (QS), have recently 
, attracted a great deal of attention in nonlinear optics, 
' where they can be easily observed experimentally in the 
T-H spatial domain (as self-supported cylindrical beams in a 
bulk medium 0|, or stripes in planar waveguides [^). 
The experimental observation of QS in the temporal do- 
main dl, and of spatiotemporal QS, or "light bullets", 
in three-dimensional samples has been recently re- 
. ported too, although it proved to be much more difficult 
to observe temporal solitons than the spatial ones. 
^ Besides optics, quadratic nonlinearities coexisting with 

'"^ linear dispersion (which is necessary to support solitons 
^ \ in balance with the nonlinearity) occur in many other 
O ' physical media, which suggests a possibility of theoreti- 
' cal prediction and experimental observation of QS in sys- 
^ r tems other than optical. In particular, it is expected that 
the experimental observation of QS in the temporal and 
spatiotemporal domains, which is fairly difficult in the 



' optical media, could be easier in other physical settings. 
In line with this idea, temporal and spatiotemporal QS 
have been recently predicted in a system of resonantly 
coupled internal waves in an inviscid stratified fluid |^ . 

An objective of the present work is to predict QS in 
dynamical lattices featuring quadratic nonlinearity. The 
resonant second harmonic generation (SHG) in nonlin- 
ear lattices, as well as more general resonant three- wave 
interactions, have been described in Refs. |6l|7tl. How- 
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ever, the scaling adopted in those papers did not imply 
existence of QS. In order to obtain them, one needs to in- 
clude into consideration the dispersion (which is by itself 
an intrinsic property of a discrete system). 

We will demonstrate that any generic quadratic non- 
linearity can easily give rise to a soliton, provided that 
the model's linear part satisfies necessary matching con- 
ditions: in the lowest-order approximation, phase and 
group velocities of the fundamental and second harmon- 
ics (to be abbreviated as FH and SH, respectively) must 
coincide at some value of the wavenumber. It turns out 
that these two conditions are difficult to satisfy in some 
standard dynamical lattices. For instance, one can check 
that the matching never takes place in monoatomic lat- 
tices, while in a diatomic lattice, whose linear spectrum 
consists of two branches, the matching occurs only in two 
cases (see, e.g., 0), when the FH wavenumber is exactly 
at the edge or center of the first Brillouin zone (BZ). 
In either case, the coincidence of the FH and SH group 
velocities is trivially provided by the fact that both are 
equal to zero. Even this trivial matching in the diatomic 
lattices can give rise to QS with a zero or small velocity. 
However, our aim here is to propose a physically realistic 
model, describing a chain of two-component dipoles in an 
anisotropic host medium, that readily meets the match- 
ing conditions in a nontrivial way (while a continuum 
limit of the model can do this only in a trivial degener- 
ate situation). We argue that the system may, in fact, 
be a representative of a wider class of models (e.g., those 
describing long molecular chains), where solitons can be 
generated by a nontrivial resonance between FH and SH. 

In section 2 we formulate the model and produce re- 
sults of its linear analysis. The nonlinear equations for 
the slowly varying amplitudes of the resonantly interact- 
ing FH and SH waves are derived in section 3. A gen- 
eralization to the two- and three-dimensional cases are 
briefly considered in section 4. In the one-dimensional 
case, the equations turn out to be essentially equivalent 
to the standard SHG equations in nonlinear optics that 
are well known to give rise to stable QS (see, e.g., Ref. 
1^ ) , while in the multidimensional case they may be more 
general than those known in optics |^ . The paper is con- 
cluded by section 4. 

II. THE MODEL AND ITS LINEAR ANALYSIS 

We consider a chain of dipoles with two components 
Xn{t) and Unit) immersed into a nonlinear anisotropic 
host medium. Physically, the dipoles may represent, e.g., 
bound electrons in a chain of impurity atoms, or in a 
string of atoms which is a part of the host medium. 

Because of the anisotropy of the system, the intersite 
coupling constants K^.y of each of two sublattices, Xn(t) 
and yn{t), may be different. As usual, we assume that 
the medium's fundamental symmetry group lacks inver- 
sion invariance, hence terms quadratic in the dynamical 
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variables may appear in the equations of motion. Thus, 
the model takes the form 



— O'xx'^n H" ^xy'^nVn H~ ^yyVn ' (2.1a) 

jjn + ^IVn - Ky (y„+i - 2y„ + yn-l) 
^ b^xx'i + bxyXnVn + byyul , (2-lb) 

where Vl^^y are the frequencies generated by the onsite 
potential, the coefficients a and b describe the onsite non- 
Hnearity, and the effective masses of the two components 
do not appear in the equations as we can divide by them 
and accordingly renormalize the other coefficients. The 
form in which these equations are written implies that 
lattice spacing is set to be equal to one. 

Note that cubic nonlinear terms, that were omitted in 
Eqs. (1) and (2) with the only purpose to render the 
equations shorter, can be readily added. It is well known 
that adding the cubic terms do not, generally, destroy QS 
1^,0. It is also relevant to note that, although we wrote 
only onsite nonlinearities in Eqs. (1) and (2), intersite 
nonlinear terms (accounting for an anharmonism of the 
intersite couplings) may also be added to the model. The 
final results presented below in sections 3 and 4 are not 
essentially altered by the intersite nonlinearities. 

Linearization of Eqs. (1) and (2) and substitution of 
XmUn ~ exp(ifcn — iojt) yields two disjoint branches of 
the dispersion relation uj{k) for the two components: 



{k)^^l,y + ^K^,y^^^^^ ik/2). (2.2) 



Note that, in principle, the coupling constants K^^y may 
be both positive and negative, corresponding, respec- 
tively, to the intersite repulsion and attraction (in the 
case of the chain of dipoles, the repulsion and attraction 
correspond, respectively, to the dipole components trans- 
verse and longitudinal with respect to the chain, so that 
at least one coupling is always positive, because at least 
one component is transverse). In any case, we assume 



that both dispersion branches (2.2) are stable, which re- 
quires |k| < il^, in the case when the coupling is negative. 

We are looking for a case in which the dispersion 
branches generated by Eqs. (1) and (2) correspond, re- 
spectively, to FH and SH. Then, the matching conditions 
necessary for the existence of QS read that, at some point 
k = ki, 



and 



ujy{2ki + Q) = 2ujx{ki) + Auj, (2.3a) 



,{ki)=Vy{2ki+Q), (2.3b) 



where Vx.y{k) = diOx.y/dk are the group velocities of the 
respective modes. Aw is a small frequency mismatch, 
I Aw I <C '^x,y, and Q is a vector of the reciprocal lattice: 
at Q = and Q — ±27r, one has normal and umklapp pro- 
cesses, respectively. The presence of Q in the matching 
conditions is a specific feature of the discrete system. 
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First we concentrate on the exact matching conditions, 
when Aw = 0. After a straightforward algebra, these 
conditions yield 

cos fci — Kx/ny , (2.4) 

hence the matching demands that {nx/ny)'^ < 1, i.e., the 
fundamental harmonic must be excited in a sublattice 
with a weaker intersite coupling. With the wavenumber 



being fixed by Eq. (2.4), another final condition gener- 
ated by the matching equations (2.3a) and ( 2.3b| ) is a 



constraint imposed on the model's parameters, 

ml -nl^ 4k-' (Ky - . (2.5) 

Of course, QS are possible not only when this constraint 
is satisfied exactly, but also in a vicinity of the corre- 
sponding surface in the system's parameter space. 

At the point where the matching conditions are exactly 
satisfied as described above, the FH frequency can be 
easily found as uji = Wj;(fci), 

Luf ^ Ql+ 2KxKy' (Ky - Kx) , (2.6) 

and the SH frequ ency is exactly 2uji 



Eqs. (2^) and (2_^) are the main results obtained from 
the analysis of the model's linear part. It is interesting to 
note that the nontrivial solution to the matching equa- 
tions is specific for the discrete model considered. Indeed, 
straightforward consideration of its continuum counter- 



part (which can be formally obtained from Eqs. (2.2) in 



the limit of infinitesimally small k) easily shows that the 



matching conditions impose not one (cf. Eq. (2.5)) but 
two constraints on the model's parameters, viz., — i^y^ 
and Q,y — 2^x- If these two conditions are met, the 
SHG resonance in the continuum model becomes degen- 
erate, holding at all the values of k. Thus, the case when 
the SHG resonance is possible in the continuum model is 
quite trivial. 

The model (1), (2) assumes no linear coupling between 
the two sublattices. However, if a weak linear coupling 
between them is present, it may be treated as a small per- 
turbation. Analysis shows (we do not display here boring 
technical details) that a weak linear coupli ng b etween x„ 



and Un only slightly changes the relations (2.4) and (2.5). 
Thus, these results are structurally stable and can be ex- 
tended to more complicated models. This, in particular, 
means that the QS solitons may also be generated by the 
SHG matching in molecular systems described by two 
linearly coupled one-dimensional lattices (see, e.g., Refs. 
II)- 



III. NONLINEAR ANALYSIS 

It is customary to seek for QS by means of multiscale 
expansions, which in our case must be implemented in 
the discrete model. To this end, we introduce a formal 
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small para mete r /i (/i <C 1), and look for a solutfon to 
Eqs. (|2.1a|) as (|2.1h| ) in the form 



MM,(ii,...;a,...)e*™* + c.c. (3.1a) 



l?Ay{h, + c.c. , (3.1b) 



where we have introduced a set of spatial and temporal 
variables tj = fi^t and S,j = ^^■'n [j ~ 1,2, ...) which are 
regarded as independent and continuous ones, A^^y are 
amplitudes which, in principle, depend on all the vari- 
ables, and c.c. stands for the complex conjugate. It is 
assumed that the carrier wave number and frequency, 
k and u j, ex actly satisfy the SHG matching conditions 
( pa|) , ( |l3b| ). 



We notice that the above substitution implies an in- 
finite lattice, which, however is not a limitation for the 
analysis: essentially the same results can be obtained for 
a finite lattice with a large numb er of sit es, see e.g. ||l3|,|7| . 

Su bstitu ting the ansatz ( 3.1a ), ( |3.1b| ) into Eqs. ( |2.l£ 



and (2.1b), and collecting all the terms of the same order 
with respect to n, at the lowest (second) order we recover 
the dispersion relation, and then at the third order we 
obtain that Ax^y depend on ti and only through R = 
^1 — vti , V being the group velocity of the modes. Finally, 
at the forth order with respect to the small parameter /i 
we arrive at a system of continuous equations. 



.dAy ^ uj'^d^A. 



^ J- XyAy^'^-'^ = (3.2b) 



dt2 2 9i?2 



where = (Pujx/dP\k=ki and ^2 = d'^^^x/dP\k=2k2+Q 
are the group velocity dispersions of the first and second 
modes, Xx = a^x/i'^i^i), Xy = bxx/i^uJi), it is assumed 
that Auj — y?6uj with 5uj ~ 0{uji), so that Alj • t = 
6uj • t2, the overbar stands for the complex conjugation, 
and dependence of A^.y o n ^2 is neglected. 

Eqs. ( |l2a| ) and ( pb| ) are identical to general equa- 



tions describing SHG in the temporal domain in disper- 
sive optical media (see, e.g., [||), 5u) playing the role of 
the standard mismatch parameter. Thus, these equations 
give rise to stable solitons of the well-known form. 



IV. MULTIDIMENSIONAL CASE 

Coming back to the interpretation of the underly- 
ing model (1), (2) as that describing interaction of 
two-component dipoles pinned at the lattice sites, it 
is straightforward to extend it to the two- and three- 
dimensional (2D and 3D) cases. Because the host 
medium is assumed anisotropic, in the most general case 
all the coupling constants may be different. Thus, in the 
2D case, for instance, the dynamical variable are Xnm and 
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Unm, the "discrete Laplacian" in Eq. (1) being replaced 
by 

where Ati'^ and Ki*"* are the coupHng constants in the 
longitudinal and transverse directions. Eq. (2) is to 
be modified in a similar way. Note that solitary waves 
in the corresponding monoatomic lattice in the strongly 
anisotropic limit, k^"^ <C Kx\ have been considered in 
Ref. 0. 

In the 2D case, the matching conditions read: 



and 



Wj,(2ki + Q) = 2cj^(ki) + Acj, (4.2) 



v,(ki) = Vj,(2ki +Q) (4.3) 



where v^.yCk) — dujx,y{'k) /dk, the bold letters standing 
for 2D vectors having longitudinal and transversal com- 
ponents; for instance, k = (fc'-'\ k^^^). 

In the 2D version of the present model, the exact 
matching conditions can be cast into a form 

(1) (t) 

1,(1) _ '^^ _ '^a; 

cos fCi 7TT" , COS fC-i 77T" . 

y y 

More sophisticated schemes with a 2D configuration of 
the FH and SH wavevectors are possible too, but we do 
not consider them here. 

The derivation of the equations for the envelope func- 
tions, outlined above for the ID case, can be easily ex- 
tended to the 2D case. As a result, at the second order 
in the formal small parameter fi defined above, the result 
is 

^ + (vVi)A.,, = 0, 

where Vj = d/d£,j and £,j is now a 2D vector, = 

From here on we conclude that the envelope 
amplitudes are functions of the traveling- wave argument 
R = ^1 — vti. Finally, at the third order in /x, we obtain 
a system 

dAx 1 V- 
i-^ H :r > 



dt2 2 ^ 

ai ,a2— l,t 



-fX.^xA^e'*'"*^ = 0, (4.4a) 



. dAy 1 ^ d'^Ay 



dt2 2 ^ "'"^'"^^i^^iai?"! 



-rXyAle^''^'^ = 0, (4.4b) 

where iOa,ax,a.2 — 9^a;a/9/c"i9fc"^ |k=ki (a = x^y and 
ai,Q;2 = (1), (t)) is the group-velocity-dispersion tensor. 
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The generalization for the 3D ca se is straigh tforw ard. 
It is necessary to stress that Eqs. ( 4.4a ) and ( 4.4b ) are 
more general than the multidimensional SHG equations 
known in nonlinear optics |9| , as they are fully anisotropic 
and contain mixed second derivatives in the variables Rj . 
Accordingly, these equations can give rise to more gen- 
eral multidimensional solitons than those known in the 
optical models However, a detailed consideration of 
these more general multidimensional solutions is beyond 
the scope of this work. 



V. CONCLUSION 



In this work, we have proposed a two-component lat- 
tice model, in its one- and multidimensional versions, 
which may satisfy the matching conditions necessary for 
the generation of solitons through the second-harmonic 
generation. The model describes an array of linearly cou- 
pled two-component dipoles in an anisotropic nonlinear 
host medium. Unlike this discrete system, its contin- 
uum counterpart meets the matching conditions only in 
a trivial degenerate situation. We have derived nonlin- 
ear evolution equations for slowly varying envelope func- 
tions of the resonantly coupled fundamental and second 
harmonics. In the one-dimensional case, these equations 
naturally prove to be tantamount to the standard system 
well-known in nonlinear optics, while in the multidimen- 
sional case the system is more general than the optical 
one, thanks to the anisotropy of the underlying lattice 
model. 

The analysis developed in this work can be extended in 
different directions. In particular, one can consider more 
general three-wave resonant interactions, which are also 
known to give rise to solitons in optics . In the context 
of condensed-matter lattice models, resonant three-wave 
interactions are, e.g., known in the model of the electron- 
phonon system of the Davydov's type , but an analysis 
shows that the full set of the matching conditions in that 
model can only be met in a limiting case of equal maxi- 
mum electron and phonon frequencies, which is not phys- 
ically realistic. A fully matched three-wave resonance is 
possible in other models, but a detailed consideration of 
this problem will be presented elsewhere. 

Note also that, in the model considered above, the cou- 
pling of the two components is purely nonlinear. Adding 
a linear coupling results in a more complicate dynamics, 
when QS are excited together with so-called accompany- 
ing modes (i.e., nonresonantly excited ones). 

Lastly, dynamics of real crystallic lattices is usually 
characterized by conspicuous dissipative losses. The sim- 
plest way to compensate the losses is by an external ac 
drive applied to the edge of a sample, the drive's fre- 
quency being equal to the FH frequency. However, de- 
tailed analysis of the loss compensation by an external 
pump should be a subject of another work. 
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